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Abstract. Using an improved version of the Hartree approximation, allowing for 
ensembles of inhomogeneous configurations, we show in a A</> 4 theory, that initially 
the system thermalises with a Bose-Einstein distribution. For later times and larger 
couplings we see deviations. 



INTRODUCTION 

In many parts of high energy physics like early universe physics and heavy ion 
collisions one wants to follow quantum field systems in time. Typical phenomena 
include reheating after inflation: after a phase of parametric resonance most energy 
is in a small part of the spectrum, a non-thermal initial condition. 

In order to study these systems and to describe phase transitions and thermal 
fluctuations between classical minima, non-perturbative approximation schemes are 
needed. Examples of high temperature approximations, that have been used suc- 
cessfully are the classical approximation [1,2] and the Hartree or gaussian approxi- 
mation [3] . The latter assumes a gaussian density matrix, such that all information 
is contained in the 1- and 2-point functions, higher point functions can be factorised 
using Wick's theorem. 



INHOMOGENEOUS HARTREE ENSEMBLE 

Instead of the commonly used Hartree approximation, which has a spatially con- 
stant mean field we used an improved version, allowing for inhomogeneous config- 
urations. We studied this approximation in a scalar ct> 4 theory in 1 + 1 dimensions: 

C = —d^dy - \^ - (1) 
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In the gaussian approximation we can expand the operator field in mode functions: 



}(x, t) = <f>(x, + E (aafa(x, t) + alf*(x, tj) , (2) 

a 

where the time-independent a a and a) a satisfy the usual commutation relations. We 
can choose them such that: 

(a ] a a a ) = n°. (3) 

For a free thermal system the n° a have a Bose-Einstein form and the f a (x,t) are 
plane waves, with suitable normalisation. 

The Heisenberg equations of motions in the gaussian approximation 2 become: 

[<9 2 -A + /i 2 + \<p 2 (x,t) + 3XE(x,t)] <f>{x,t) = 0, (4) 
\d 2 t - A + /i 2 + 3\<p 2 (x, t) + 3AS(x, t)]f a (x, t)= 0, (5) 



E(M)=E |/«(x,0| 2 (2n°+l). (6) 

We mention that there exists an effective Hamiltonian which leads to the above 
equations, suggesting that in the end the system will equilibrate classically accord- 
ing to this H c k. But we hope that this has a timescale much larger than that of all 
interesting processes. 

Initial conditions We need to specify the initial fields and their conjugate 
momenta, which together with the n° a specifies an initial gaussian density matrix. 
By averaging over different runs with different initial conditions we can build non- 
gaussian density matrices, allowing for much more general initial conditions. Only 
the time evolution is evaluated with gaussian density matrices. 

We will take the mean field in its zero-temperature minimum, with all energy in 
a few of its momentum modes, like after parametric resonance: 

Jmax {X \ 

(f)(x,0)=v and <fr(x, 0) = A ^ cos ( 2irj— - ipj ) , (7) 

where the ^-'s are random phases and 27r? max /L is of the order m. We will average 
over several runs, typically 10-20. 

For simplicity we will take the quantum modes equal to the zero-temperature 
vacuum form: 

2 ) The classical approximation can be obtained by the substitution 3 — > in Eq. (4) while 
omitting Eq. (5). The large n approximation leads to 3 — > 1 in Eq. (4) and Eq. (5). We chose not 
to use large n to avoid problems with would-be Goldstone bosons in 1 + 1 dimensions. 
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FIGURE 1. Left: time dependence of the average mean field. Right: the energy densities for 
the same run: total (straight line at 1), mean field (starting at total energy), and modes. 



with ujk = \Jm 2 + k 2 , m chosen as the zero-temperature mass and L the size of the 
system. To put the system on a computer we discretise it on a space-time lattice. 
In the limit of (spatial) lattice spacing a to zero, the mode sum E occurring in the 
equations of motion is logarithmically diverging, but this can be fixed by a simple 
mass renormalisation: fi 2 — > ji 2 cn = jj? — |£ log(Aa 2 ). 

Thermalisation If A is small and the system is not too far from equilibrium 
the interaction of the mode functions with the typically inhomogeneous mean field 
can be viewed as scattering of particles via the mean field. Because of this scat- 
tering they will hopefully thermalise according to a thermal - i.e. Bose-Einstein - 
distribution. 

To check this, a definition of particle number is needed. For a free thermal system 
the 2-point functions read: 
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(9) 



We want to define an rik and Uk in a similar way as Eq. 9. This can be done 
by coarsening the measured 2-point functions over the initial conditions, a time 
interval, and space. 



NUMERICAL RESULTS 

In Fig. 1 we plotted for a typical run (256 lattice points, size Lm = 32, \/m 2 = 
1/12) the spatial averages of the mean field and the energy densities as a function 
of time. We see that the energy, initially completely in the mean field, is going 
to the modes with a timescale of the order 100 in mass units. The total energy 
is conserved up to a few percent. The energy staying in the mean field - around 
20% - is much more than expected from classical equipartition: around 1/256 
for 256 mode functions. Fig. 2 shows a plot of the calculated particle numbers 
obtained after averaging over whole space, approximately 5 oscillation periods and 




10 configurations. On the left the situation is plotted at shorter times, on the right 
at larger times. We see that starting at low momenta, a Bose- Einstein distribution 
is emerging. Only the zero mode is showing deviations. 

Runs at twice this energy density show very similar results, but with faster 
thermalisation. At large times deviations from the straight line BE's start to show 
up, which we interpret as signs of classical equipartition. For strong coupling or 
after a long times we see the emergence of an offset, which might be a chemical 
potential following from approximate particle conservation. 

We finally mention checks with different kinds of initial conditions (not presented 
here), which also show thermalisation according to a Bose- Einstein distribution. 



CONCLUSION 

Using our Hartree ensemble approximation with inhomogeneous mean fields, we 
find thermalisation according to a Bose-Einstein distribution, starting in the low 
momentum modes. The time-scales heavily depend on the size of the coupling and 
the energy density. 

Signs of classical equilibration occur only at timescales much larger than that for 
the onset of the Bose-Einstein thermalisation, so even if they will prevent the BE 
thermalisation to reach arbitrary high values of u>, they do not influence processes 
with an intermediate timescale. 

An extension to higher dimensions is wanted but numerically demanding. A 
possible way out is a reduction of the number of mode functions. 

A more elaborate discussion of these simulations is in preparation. 
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